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Recent radio VLBI observations of the complex gravitationally lensed system 
PKS1830-211 have thrown up some questions with regard to the processes occur- 
ing at the heart of the blazar source at a redshift of 2.51, which is viewed almost 
straight down the jet axis. This work links, by a model of a helical jet tracked by bal- 
listically ejected plasmons from a precessing nozzle, observations on the scale of tens 
of microarcseconds to those on the scale of milliarcseconds. An observed jet precession 
period of 1.08 years is inferred from the model, translating to an intrinsic period of 
30.8 years for a source at redshift Zg = 2.51 and an assumed jet bulk velocity P of 
0.99c. This fits well with the picture of the active galactic nucleus hosting a binary 
black hole system at its centre, with the jet emitted by one member of the system, 
and the precession as being due to its orbital motion around its companion. 
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INTRODUCTION 



The much-studied gravitationally lensed system, PKS1830- 
211, was identified as such by Rao & Subrahmanyam (1988), 
and is one of the strongest radio sources in the sky at cen- 
timetric wavelengths. On arcsecond scale, it consists of two 
prominent radio core-knot structures, lensed images of each 
other, each core-knot axis being roughly perpendicular to 
the line of separation of the two cores (the two cores be- 
ing separated by about .98 and aligned roughly northeast 
(NE) and southwest (SW) in the plane of the sky). The 
two core-knot groups (hereafter NE and SW groups) are 
point-inversion symmetric with respect to the centre of the 
system, and show somewhat weaker diffuse emission (Rao 
& Subrahmanyam 1988, Subrahmanyam et al. 1990); in the 
radio L-band, the diffuse emission is extended enough that 
its images actually close round, linking the image groups 
and making a pseudo-Einstein Ring system (Jauncey et al. 
1991), with nearly half of the total flux density of about 12 
Jy coming from the cores. The cores themselves have flat ra- 
dio spectra (Rao & Subrahmanyan 1988). On the scale of 
milliarcseconds in the radio at 43 GHz, a weaker, inner knot 
accompanies each of the bright cores, on the same side of 
each core as the arcsecond-scale knot (which is not detected 
at this resolution), and separated from it by ~ 1 mas (Gar- 



rett et al 1997, Jones et al., Guirado et al. 1999, Jin et al. 
2003). 

Jin et al. (2003) infer that there are time-dependent po- 
sitional variations in the radio centroids of emission in the 
bright cores at 43 GHz, resulting in changes in the separation 
of the two core components on the scale of lO's to lOO's fias 
(Figure between successive observational epochs. Signifi- 
cant radio variability has been reported (Lovell et al. 1996, 
van Ommen et al. 1995, Hagiwara et al. 1996 and Lovell et 
al. 1998); in fact, Lovell et al.(1998) report a time delay be- 
tween the arrival time at the observer of variations in one 
core image relative to (correlated) variations in the other, of 
26I5 days. However, Romero et al.(1997) report this system 
to be non- variable, based on data taken during 1995-1996 in 
the radio L-band. Lovell et al's (1998) analysis was based on 
8.6 GHz data between 1996 and 1997, during which period 
dramatic variations occured. 

The source has been identified in EGRET observations 
(Mattox et al., 1997). 

In the X-ray, ROSAT obseravtions have revealed ab- 
sorption against the compact images, in excess of what may 
be reasonably attributed to our Galaxy, implicating a fore- 
ground (lensing) galaxy (Mathur & Nair 1997). ASCA obser- 
vatons (Oshima et al. 2001) suggest significant X-ray vari- 
ability in the system, so that the ratio of the NE image's 
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emission to that of the SW image was a factor of 7 (as com- 
pared with a typical factor of 1 to 1.8 in the radio). Given 
that the source properties resemble those of a blazar, this is 
very likely a variation intrinsic to the blazar rather than a 
lens-induced one, though observations of correlated spectral 
hardness changes would be necessary to confirm that such a 
phenomenon indeed occurs in this system. 

The system that acts as the gravitational lens, now 
known to be fairly complex, is believed to be in the main a 
gas-rich massive spiral galaxy, with a measured redshift of 
zi = 0.89, through molecular absorption line studies (Wik- 
lind & Combes 1996, Carilli et al. 1998, Chengalur et al. 
1999). In the optical, observations are hindered by the loca- 
tion of the system (in Galactic coordinates, I = 12°. 2, b — 
—5°. 7; here optical extinction is experienced of order 2.7 
magnitudes, cf. Subrahmanyan et al 1990); nonetheless Frye 
et al. (1999) (I and K band observations), Courbin et al. 
(1998, 2002) (V, I, J, H & K bands), Lehar et al.(2000) 
(H and I bands) and Winn et al. (2002) (1 and V bands) 
have presented a wealth of optical and near-IR data, which, 
suitably image processed, reveal a spiral galaxy lens seen 
almost face on, with possibly a secondary lens nearby, and 
perhaps some sign of a lower redshift spiral galaxy about 
2 .5 southwards of the system. Lidman et al. (1999) deter- 
mine through infrared spectroscopy a redshift for the source 
{zs = 2.51). HI observations (Lovell et al. 1996) have earlier 
revealed some absorption at redshift Z2 = 0.19. The presence 
in the optical and IR bands of a foreground low redshift spi- 
ral galaxy is highly suggestive of the existence of a secondary 
lens. 

Gravitational lens modeling of the system that has at- 
tempted to reproduce a range of image properties (Subrah- 
manyan et al. 1990, Kochanek & Narayan 1992, Nair et al. 
1993), has shown the system to be consistent with an iso- 
lated elliptical lens of mass of the order of that of a spiral 
galaxy. However, in view of the apparently complex opti- 
cal/IR structure subsequently observed along the foreground 
and upto the lens redshift, a detailed remodeling is in order. 

The present work seeks to provide an explanation for 
the puzzling observations of Garrett et al.(1998) and Jin et 
al.(2003), wherein it is shown that structural and temporal 
variations occur in and around the radio core images in this 
system in VLBA observations at 43 GHz. These observations 
probe the system at the highest resolution to date, of consid- 
erable interest because of the combination of a high redshift 
source {zs = 2.5) and enhanced resolution due to magnifi- 
cation by lensing. For the purpose of the present study, a 
full reworking of the lens model is actually inessential; in- 
stead, certain lens properties applicable to this system will 
be either invoked or derived as required. 



2 SOME PROPERTIES OF THE LENS 

Considering the optical and infrared observations of possi- 
ble multiple lensing bodies along the line of sight to the 
system (as described in the previous section) , only the most 
elementary assumptions regarding the lens system will be 
made. At the location of each image in the plane of the sky, 
it is assumed that it is possible to define a 2 x 2 source-to- 
image coordinate transformation matrix that is symmetric. 
This is valid provided one or the other of two circumstances 



occurs: (a) the dominant lensing agents are roughly at the 
same distance from the observer, or from the source, so that 
they can be said to inhabit the same lens plane ( cf. Bland- 
ford & Narayan 1986), or (b) save for the main deflectors in 
the same lens plane, lenses at other redshifts provide only 
a shear perturbation (however strong or weak) to the main 
lens's action. Then the lensing equations reduce to one iden- 
tical with single-plane lensing, and source-to-image transfor- 
mation matrices are symmetric as a result (Kovner 1987). 
The latter assumption is adopted in view of the possible ex- 
istence of a lens ai z = 0.19, which is seen 2. 5 south, in 
the plane of the sky, of the lensed images in the system. In 
fact, since it has proven fairly successful to model this sys- 
tem with a single elliptical lens (Subrahmanyan et al. 1990, 
Kochanek & Narayan 1992, Nair et al. 1993), it is possible 
that the main lens, a spiral galaxy at zi = 0.89 which is 
viewed almost face-on, is aided by shear roughly along the 
north-south direction by the foreground lens, mimicking the 
behaviour of a single elliptical lens with such an orientation 
(note that this is close to the direction of the major axis 
of the elliptical lens model in Nair et al. 1993). That the 
effect of the foreground (z2 — 0.19) lens can be treated as a 
perturbation on the main lens' action derives some support 
from the fact that lens models place the effective lens centre 
between the NE and SW groups of images, and somewhat 
closer to the negative-parity SW group of images (the parity 
of which is demonstrated by the fact that temporal varia- 
tions in its core's radio flux density lag behind those in the 
NE core; cf. Lovell et al. 1998). This would be as expected 
for a single lensing galaxy, suggesting that the actual lens is 
not very far from one. 



3 A LENSED HELICAL JET 

3.1 Inferred source evolution from observations 
on the scale of tens of microarcseconds 

Figure0(from Fig. 3 of Jin et al. 2003), plots, for 8 successive 
epochs of VLBA observations at 43 GHz (0."3 mas FWHM 
beam), changes in declination versus changes in right ascen- 
sion (RA) of the position vector from the centroid of emis- 
sion of the NE image core to that of the SW image core, 
all relative to the first epoch (trcf) of observation. Formally, 
the error in estimating these changes is of the order of 3 fias 
as both the core images appear in the same field of view in 
each observation (Jin et al. 2003). The plot is displayed here 
with straight-line fits to two remarkably linear tracks that 
can be identified in the plot. In this section, an explanation 
for this phenomenon will be sought. 

Anticipating a physical model to be discussed later in 
this section, it is assumed that the source, a blazar, emits 
plasmons, one by one, from its central engine, and that the 
observed (but essentially unresolved) core radio emission is 
dominated for a while by that of the most recently emitted 
plasmon, which propagates away from the central engine 
with constant velocity before fading from view (its spectral 
evolution, synchrotron self-absorbed, being as described in, 
for example, van der Laan 1966). Any changes in the lo- 
cation of the radio emission of the central engine itself are 
assumed to be undetectable as a function of time (< 3 ^sls). 
(The actual position of the central engine then drops out of 
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Figure 1. 

Changes in the centroids of radio emission in the NE and SW 
cores; their separations as a function of epoch relative to Epoch 
1 in the figure (from Fig. 3 of Jin et al. 2003). Superposed are 
straight Une fits to the two Unear runs apparent in the data, which 
yield slopes of (for epochs 2-3-4) 0.85 it 0.07 and (for epochs 5-6- 
7-8), 4.66 ±0.48. 

the following computations, because we consider changes in 
position relative to a reference epoch). 

The radio core of the source and its associated features 
are lensed into two images, denoted by the subscripts A (the 
NE group) and B (the SW group) in what follows. Changes 
occuring in the unlensed source shall be referred to as being 
in the source plane, and the corresponding effects in the im- 
ages shall be referred to as occuring in the image plane. For 
the unlensed source, let the angular distance of an ejected 
plasmon from the central engine at any given epoch of ob- 
servation, tk, be denoted by A$[tfc], its components along 
the directions of increasing RA and declination being repre- 
sented by A<[>a and A<I>5 respectively. 

Now, local to a point in a gravitationally lensed image, 
positional changes are mapped through a linear transforma- 
tion of the underlying changes in the source. That is, for 
small changes A<1? in the source, the corresponding changes 
in images A and B are given by: 



A4.0 



A$o 



A$5,a) = 
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/ A$c 


V 


04 J 




bi 




( A*, 


q 




I A$,- 



(1) 
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Two comments are in order at this point: (a) It has been 
assumed that the source-to-image transformation matrices 
are symmetric, under the circumstances discussed in Sec- 
tion|21and (b) , when considering observed changes in the im- 
ages, it is necessary to account for the fact that events in the 
SW image are behind those in the NE image, by an amount 
equal to the time delay between the images, r = 26^4 days 
(Lovell et al. 1998). Thus the observed changes in the two 
images map back to two different epochs in the source itself 
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Figure 2. 

Underlying model for the processing jet. The jet emission vector 
makes an angle <j) (jet half-opening angle) with its axis, which is 
also the line of sight by assumption (left). From time to time, 
plasmons are ejected, and each propagates with constant velocity 
along the surface of the cone defined by the motion of the emission 
vector. (Right): The view along the line of sight to the blazar, 
prior to lensing. 

(epochs referring to events in the source will hereinafter be 
called source epochs, as distinct from observer epochs). 

Formally, it is possible to define a relative or image-to- 
image transformation matrix (in this case, the SW image- 
to-NE image, or 'SW2NE', matrix wiU be defined): 



(A<l>o 



A$ 
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( A$c,B \ 


[ T3 




I A$i,B J 



(3) 



where rj = {bibi-q^), Ti = {aibi-pq) /ri, T2 = {bip-aiq)/r], 
T3 = {pbi — qa4)/ri, and T4 = (0461 — pq)/v- When evalu- 
ating this matrix from the observations, which here include 
evolving components, care must be taken that appropriate 
corrections are introduced for the differing source epochs 
(or alternatively, fight arrival time delays) corresponding to 
each of the images. 

Adopting the Hjellming & Johnston (1981) model for 
SS433 (or the Gower et al. 1982 model, as applied to AGNs) 
for a precessing helical jet, we assume that a plasmon in the 
jet is ejected from the central engine of the source with con- 
stant velocity v along the surface of a cone of half-opening 
angle (f>, the symmetry axis of which is the axis of the jet 
(c/. FigHJ . New plasmons are ejected ballistically from time 
to time. The plasmon ejection velocity vector precesses with 
angular velocity Q about the cone's surface, making a con- 
stant angle (j) with the jet axis, so that each subsequent 
plasmon is ejected at a different phase of the jet's preces- 
sion. Once ejected with a given velocity v, a plasmon con- 
tinues to propagate away from the central engine with that 
velocity (or it fades from view before it appears to deceler- 
ate or accelerate). Each plasmon is assumed to consist of an 
isotropically expanding cloud of synchrotron-emitting elec- 
trons, so that the centroid of emission is not affected by its 
evolution. As mentioned earlier, it is assumed that only one 
plasmon at a time dominates the unresolved radio emission 
from the radio core. It is also assumed that the fine of sight 
of the observer coincides with the jet precession axis, for 
simplicity of calculations. 

Let a plasmon be ejected from the central engine of 
the source. The observable proper motion of the unlensed 
plasmon from the time of its ejection, t^, to the epoch, k, 
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of observation, tk, would be given by its components in the 
plane of the sky: 



Va {tk — tcj ) 

d{l~-v^/c) 

Vs{tk ~ tcj) 



(4) 



(5) 



d (1 - v^/c) 

Here d is the angular diameter distance to the source, and c 
is the velocity of light. The source time interval {tk — t^) is 
corrected for the motion of the plasmon at relativistic speeds 
along the line of sight during its propagation. (To convert 
a source time interval itself to the observed time interval, 

((5t)obscrvor = ((5t)source/ ( 1 - Va: / c) / {1 + Zs), where Zs IS the 

redshift of the source; here one corrects for cosmological ef- 
fects as well. This term is included in the angular diameter 
distance in Eqns.^l&Eland therefore need not be specified 
explicitly) . 

At epoch tk, we observe in image A: 

A'l'cA [tk] = aiA^^+pA^f 

A^s,A [tk] = P A^t + aiA^f 

where the superscript 'A' is used to indicate that an ap- 
propriate correction for the light travel time along the path 
through image A must be applied. For image B: 

A-l-.^.B [tft] = 6iA$^-hq A*? 



A$i 



[tk] 



q A$^ + 64A$f 



with a similar use of superscript as before. The convention 
that will be used in the following is that for events in image 
A, {tk — tcj) is written as {tk — A^cj), and for image B, as 
{tk — Btej) = {tk — Atcj — t), explicitly writing the light 
travel time difference r between the paths through images 
A and B into the expressions. 

At this point contact is possible with the observations 
in Figure What is observed in each image is a change 
A<I> from a location 6 = {6a, 9s), the position of the lensed 
central engine. It is assumed that 9 for each image remains 
fixed through the different epochs of observation; since the 
changes in FigureQare plotted relative to a reference epoch, 
the location 9 in each image will drop out of all expressions 
thereafter, hence it will be ignored. 

With this simplification, and making the appropriate 
substitutions, the following quantities are constructed: 



5$c [tk] = A$c,A [tk] - A<1>„,B [tk] 

= {(ai - bi)v'a + {p- q)v's}{tk - Atcj) + 

{biVa+qv's)T (6) 

S^S [tk] = AcE>«,A [tk] - A$5,B [tk] 

= {{p - q)v'a + (a4 - b4)vs} {tk - Atej) + 

{q v'a + b4Vs)T (7) 



In the above, the prime on quantities Va and vs is used to 
denote the corresponding velocity components divided by 
d (1 — Vx/c). Also note the implicit assumption that at the 
epoch of observation, tk, the plasmon responsible for the 



observations is present in both images (i.e., tk > A^oj + t). 
As long as only one plasmon at a time dominates the source 
emission on this scale of observation, and it expands isotrop- 
ically (so that the centroid of radio emission is not affected 
by its evolution), it is possible to neglect the effects of evo- 
lution of the plasmon itself over the timescale of observation 
and deal simply with positional changes of its centroid rela- 
tive to the core. The presence of a second evolving plasmon 
plasmon coexisting with the first would considerably com- 
plicate this simple picture! 

The changes plotted in Figure are relative to a refer- 
ence epoch trof, so we construct new quantities, which are 
now the abcissas and ordinates of the plot in Figure Q: 

5{5^o.) [tk] = S^c [tk] - 5^c. [tret] 

= {ai - bi) v'a + {p - q) v's (8) 
S{S$s)[tk] = S<i>s[tk]-5<i>s[trc!] 

= {P- q) + (14 - 64) v's (9) 

(10) 

The slope m of a particular linear track in Figure Q is 
given by (dropping the primes on and vs as unnecessary 
here) : 

5{5$a) [tk]-S{5$a) [tj] 



5(5$,) [t,] - 5{S<i>s) [tj] 

{a-i - bi) Vg + {p - q ) vs 
{p — q ) Va + {a4 - 64) vs 



(11) 
(12) 

From the above expressions it is now apparent that the 
change in slope of the tracks between the run of epochs 2, 3 
& 4 (slope mi = 0.85 ± 0.07) and that of epochs 5, 6, 7 & 8 
(slope 7712 — 4.66 ±0.48) is effected by a change in the quan- 
tities Va and vs , projections of the ejection velocity vector v 
of a plasmon, in the plane of the sky. In the context of our 
model, this corresponds to the dominance of a new plasmon 
between epochs 4 & 5, which has been ejected at a phase 
in the precession of the jet which is distinct from that of 
the previously ejected plasmon. An obvious question at this 
point is, what is the actual event that causes a discontinu- 
ity between epochs 1 & 2, and again between epochs 4 & 5? 
The appearance of a new plasmon in the NE image precedes 
that in the SW image by 26 days, or nearly two successive 
epochs in our 43 GHz observations, so if the discontinu- 
ity were to be caused by its appearance in the NE image, 
there should be another such discontinuity two epochs later, 
caused by its manifestation in the SW image. The fact that 
both linear tracks are more than two epochs long suggests 
that this is not the case. The reason is somewhat subtle and 
will be discussed in Section f3.6l after deriving some quanti- 
ties in connection with the source-to-image transformation 
matrices, but for the moment we identify the discontinuities 
in these linear tracks as corresponding to the appearance 
of new plasmons in the SW image (so that, by the time a 
discontinuity occurs, the same source plasmon is imaged in 
both the NE and the SW regions). 

As an aside, notice that for equal intervals between 
epochs, the changes 55^ should be the same, along a single 
linear track. That this is only approximately true for each 
of the two tracks in Figure is an indication that emission 
near the central core is variable at some level, pulling in or 
pushing out the centroid of radio emission from time to time. 
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This results in an apparent change in speed of the emergent 
plasmon, with no significant change in direction. 

For later use, the above expression is recast in a slightly 
different form: 



- ai) 



p{l — mai) — (f — mbi) 



(13) 



where p = p/q, di = ai/p, 04 = a4/p, bi = bi/q and 64 = 
64/g. For a given value of slope m, once the lensing— related 
quantities p, oi, 0,4, &i & 64 are determined, this expres- 
sion gives the phase of precession of the jet (modulo nir 
radians, where n is a natural number), relative to a suitably 
chosen reference direction, at the time of ejection of the rel- 
evant plasmon. Thus, the difference of the phases between 
the events slopes mi and 7712, over the known difference in 
epoch between the observations, yields the angular velocity 
or rate of precession of the jet's source, Q — 



(lp2 — Ipl 



(t2 -tl) 



where: 



Ipl = —Tan 



Tp2 = —Tan 



-1 fvs_ 
-1 fvs_ 



mi 



m2\ 



(14) 



(15) 



(16) 



and ti and t2 are just the dates of the first epoch of observa- 
tion of each event (ti is epoch 2, 2 Feb 1997, and t2 is epoch 
5, 21 Mar 1997). The uncertainty in the commencement of 
each event is regarded as being uniformly distributed, be- 
tween epochs 1 and 2 for ti , and between epochs 4 and 5 for 
t2 (see Fig. 0. 



3.2 The milliarcsecond scale behaviour 

At this point, it is necessary to review some puzzling VLBA 
observations of PKS1830-211 at 43 GHz which preceded 
the eight-epoch track of observations of Jin et al. (2003), 
namely, those of Garrett et al.(1998). These observations 
revealed dramatic changes in the milliarcsecond scale struc- 
ture of the two image groups between epochs of observation 
30 May 1996 and 14 July 1996. The former epoch of obser- 
vation showed core-knot structures in both the NE and SW 
images, as for each of the epochs of observation in Jin et 
al.(2003), but the latter epoch, quite suddenly and surpris- 
ingly, exhibited rich structure in each of the image groups, 
with as many as six features in the NE image group be- 
ing potentially identifiable with corresponding features in 
the SW image group (including the basic core-knot struc- 
ture, cf. Fig.|2|l. By the time of the beginning of Jin et al.'s 
(2003) observations (19 January 1997), the system was back 
to its quiescent phase, with only a simple core-knot struc- 
ture being imaged. At some point, on a timescale of several 
weeks to half a year, a dramatic outburst of plasmons oc- 
cured, which then evolved and faded from view, and did not 
recur over the three months that followed, during which Jin 
et al.'s (2003) monitoring observations of the system were 
in progress. 

It is expected that at least some of the plasmons in the 
NE image will have corresponding entities in the SW image. 
Obviously related are features Al and BI, the core images, 
and A4 and B4, which are seen in both epochs of observation 
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Figure 3. 

(Top) : The observations of 30 May 1997 of the NE (left) and SW 
(right) VLBA images at 43 GHz and 0.33 mas circular restor- 
ing beam, showing a simple core-knot structure in each image. 
(Bottom) : The same as above, but for the next epoch of observa- 
tions, 14 July 1996, exhibiting rich structure in both images. (Both 
maps from Garrett et al., 1998.) 

in FigureOland appear to be static 'knots'. Of the short-lived 
features observed in the maps of 14 July 1996, A3 and B3 
are obviously not related to each other, from considerations 
of image parity. However, A2 and B2, A5 and B5, and A6 
and B6 are easily seen to be corresponding images, once the 
effects of the time delay between the NE and SW images 
(recalling that phenomena in the NE group of images occur 
about 26 days earlier than they show up in the SW image 
group), and the distorting effects of lensing in the two image 
groups are accounted for. It is conjectured at this point that 
feature A3 has not as yet manifest a counterpart in the SW 
group of images (it is younger than 26 days) , and feature B3 
is so aged that its counterpart in the NE group has already 
evolved below the map detection limits. 



3.3 Deriving the relative (SW2NE) image 
transformation matrix 

Assuming that the source-to-image linear transformation 
matrices remain locally valid across the entire NE or SW 
image group (on scales of order 1 mas), and that each plas- 
mon in the source is moving ballistically away from the core, 
on a straight line trajectory, it proves feasible to derive a rel- 
ative (image-to-image) 2x2 transformation matrix from the 
observations of 14 July 1996 shown in Figure |H] 

Recasting the relative image transformation matrbc in 
the notation of EQuation ll3l 



S,A ) 





T2\ 


/ A$„.b\ 






I A<1>4,B ) 



(17) 
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with Ti = K(ai64 - 1), T2 = K(6i - di), T-j, = ^(64 - 04) and 
74 = K(a4foi - 1), and N = p/(bJn - 1). 

For a given feature i, its separation from the core posi- 
tion wiU be denoted by the vector (A^a.Aj, ^'^'iS.Ai)- From 
the static features, the cores Al and Bl, and knots A4 and 
B4, one derives the core-knot separation vector in each im- 
age group, yielding two constraints (positional coordinates) 
on the relative transformation matrix (from Eqn. 1171 : 



TiA$c«.B4 + r2A$,-,B, - A$„,A4 = 
T3A-1-<,.B4 + r4A$5,B4 - A$i,A4 = 



(18) 
(19) 



The pairs of relationships between A2 and B2, A5 and 
B5, and A6 and B6 each yield a single constraint, from 
the shared direction of motion of the plasmon in the source 
(these are evolving features and their actual locations in the 
NE and SW images refer back to different source times due 
to the relative time delay). That is, given the location of, say, 
plcismon B2 in the the SW group relative to its core image 
Bl, the relative transformation matrix must be such that it 
reproduces the correct direction for plasmon A2 relative to 
its core Al in the NE group of images. 
For features A2 and B2, A5 and B5, and A6 and B6, the 
single constraint equation from each pair reads (with i = 
2,5,6): 



A$, 



S,Ai 



A$„,Ai ~ TiA-l>„,B, -f TaA^i.B. 
which may be rewritten as: 



(20) 



(21) 



TlA$c«,BiA$5,Ai -HT2A$5,AiA$5,Bi- 
r3A$c,A,A-l><,,B. - T4A<l>,,A.A<l>i,B, = 

Thus the observations yield, through Equations 1181 IHH 
and 1211 5 constraint equations for the 4 unknown relative 
transformation matrix elements, Tj, j = 1 to 4. Note that 
the unknown matrix elements and the observations deter- 
mining them are implicitly related, non-linear in the obser- 
vational quantities, and of the form: F(Xa,La) = 0, where 
Xa is the vector of 4 unknowns (matrix elements) and La is 
the vector of 16 observational quantities. F itself consists of 
5 nonlinear functions. The unknowns in vector X are solved 
for via a mixed adjustment model (see, for example, A. Leick, 
1995), the solution of which is described in the Appendix. 

The solution for the elements of the SW2NE transfor- 
mation matrix are: 



-1.279 0.871 
0.973 0.245 



(22) 



Its determinant is: —1.161. This is the NE/SW image flux 
density ratio, k — TiTi—T2T-i. The covariance (error) matrix 
is given by: 



/ 0.071 0.031 0.106 0.053 

0.031 0.017 0.054 0.026 

0.106 0.054 0.342 0.165 

0.053 0.026 0.165 0.081 



(23) 



From the expressions immediately following Equation II 71 it 
is possible to express the quantities ai, 04, 61 and 64 in terms 
of Ti, T2, T3 and T4 and p. 



a~i = (Ti - k/p)/T3 
04 = (Ti - k/p)/T2 
bi^ip- Ti)/Ti 
hi = {p-Ti)/T2 



(24) 
(25) 
(26) 
(27) 



The evaluation of these quantities depends on the value of 
p, for which straightforward algebraic methods of determi- 
nation exist, but which were found to be sensitive to numer- 
ical errors of estimate. Hence a different approach has been 
taken, which will be discussed in Section 1^31 



3.4 The directions of ejection for the 
miniarcsecond-scale features 

It is possible, under the explicit assumption that the source 
velocity for each evolving plasmon remains constant from 
the time of generation at least upto the epoch represented 
by the observations of 14 July 1996 in Figure |3 to relate 
the distance travelled by a given plasmon as seen in the NE 
image group of that observation, with the distance that is 
travelled by its corresponding feature, transformed from the 
SW group of images via the relative image transformation 
matrix Tsw2ne and corrected for the lag in time of 26 days 
for the SW group. 

That is, the squared distances are (for i = 2, 5, 6): 



D 



A* 



+ A*: 



i,Aii 



(28) 



in the NE group of images, and transformed from the SW 
group to the NE group, are: 

D|w-2iVi5= (TiA$^,B, +r2A$i,B.)'+ (29) 
(T3A$<,,B. -f r4A$i,B,)2 

In the NE group, the velocity of motion is Djv_B/(ffe — tej), 
and, transformed from the SW group, the corresponding ve- 
locity of motion is Dsw2ne /{tk — tej — t), where tk is the 
epoch of observation and r is the time delay between the im- 
age groups (assumed to be constant over the span of about 
a milliarcsecond from the core). If Q represents the ratio 
^ne/Gsw2nb, then: 

tQ 



(30) 



with Q > 1. 

In order to estimate how well the mixed adjustment 
model-derived relative transformation matrix actually re- 
produces the position angles of the features in the NE 
group of images when it is applied to data from the SW 
group, the quantities A"I>i,Ai/A"I>Q,,Ai (for the original posi- 
tion angle of feature i in the NE group), and (r3A<l>Q,,Bi + 
T4A$5,Bi)/(TiA$a,Bi -|-r2A$a,Bi) (for the position angle 
of the transformed data relating to feature i from the SW 
image) are calculated. For feature A2 (B2), the original po- 
sition angle is —48° versus a transformed value of —30° but 
Q is < 1 (pointing to some inconsistency for this feature); 
for feature A5 (B5), the original position angle is 104° as 
compared with a transformed value of 113°, and Q is 1.272; 
from Equation 1301 a value of 125.2 days is obtained as the 
'age' of the feature at the time of observation in the NE 
image on 14 July 1996. For feature A6 (B6), the original 
position angle is 159°, and the transformed position angle 
is 158°; Q is 6.255 and the 'age' of the feature as seen in 
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the NE image is 30.9 days. (Owing to the existence of signif- 
icant cross- correlations between quantities derived through 
the mixed adjustment model - c/. Equation 1231 - formal er- 
rors on quantities derived will be given only for the final 
result). A cross-check on the mapping of the static knot re- 
veals a perfect match for position. The position angle is 99° 
in the original; transformed, it is 99°, and as to the distance 
from the core, it is 0.99 mas in the original and 0.99 mas in 
the transformed version. Thus the transformation matrix is 
heavily influenced by feature A6(B6) and the the location 
of the static knot feature. 

Of course, in order to derive the rate of precession of 
the jet from these numbers, it is necessary to consider the 
mappings in the source plane of these features, rather than 
those in the NE image. For this, it is necessary to calculate 
the velocities of the various evolving features in the source 
plane. The elements of the source-to-image transformation 
matrix in Equation^can be rewritten, dividing each element 
by p. Then its inverse reads: 

Ts2JVi5 =-P~^(«i"4 - 1)"^ ^ (31) 

Similarly, for Equation |2l the matrix can be recast and in- 
verted to yield: 

T^Slsw = a"'(6i&4 - 1)-' ^ (32) 

Taking the time-derivative (denoted by a superscribed dot) 
of the angular separations involved, so that one is now deal- 



ing with velocities. 










Ac&o, — p~^[dia4, — 


1)- 




~ A$i,A) 


(33) 




1)- 


(biA'l'a.fl - 


- A'$5,b) 


(34) 


A$i = p~^{d-iaA - 


1)- 


(aiA$5,A - 


A-I-ca) 


(35) 


= q-'(6l64- 


l)-^l 




A$<,,s) 


(36) 



Now, from the inverse of the image-to-image transformation 
matrix, 

A<l>,,s = fc-'(r4A-I'„,A - T2A$i,A) (37) 

A$i,fl = fc~'(riA$i,A - TsA-I-ca) (38) 

so that, substituting for di and £4, or b\ and 64 from Equa- 
tions]^ to |57| and with a little algebra, one obtains: 

A$^ ^ T2 {PA^S,^ - A$j,A) ^ VS_ ^gg^ 

A$^ 73 (pA$„,B - A$c,a) 

This expression yields (as in the microarcsecond-scale anal- 
ysis discussed earlier) the phase (modulo n-n) of ejection of 
the relevant plasmon, with a suitable choice of reference di- 
rection. The quantities A$a,B and A<I>i,B are obtained from 
A$c«,B and A$a,B by dividing by {tk — tej — r), and A$a,A 
and A$5,A are obtained by dividing A$a,A and A"l>a,A by 
(tk — tej). To estimate the rate of precession of the jet from 
expression 1391 the phases of ejection, tpj, are calculated for 
plasmons j (with j = 5, 6) and the observed time that has 
lapsed between the ejection of plasmon 5 and plasmon 6 is 
calculated as (tfe — tej.e) — {tk — tej,5), where tk, the shared 
epoch of observation, drops out of the estimate. The rest 
of the calculation follows along the lines described for the 
microarcsecond scale analysis. 



3.5 Evaluating the period of precession of the jet 
from milliarcsecond and microarcsecond scale 
information 

Expressions 1131 and 1391 are, respectively, the phases of ejec- 
tion of various plasmons on the microarcsecond and milliarc- 
second scales, as derived from observables. Note that both 
expressions require the evaluation of p. In order to do this, it 
will be assumed that the period of precession of the relativis- 
tic jet remains the same when evaluated from data on the 
scale of tens of microarcseconds, as from data on the scale of 
milliarseconds (for a jet opening half-angle of, say, 0°.5 and 
assuming a typical linear magnification by the gravitational 
lens by a factor of 3, this difference of scale corresponds to 
source-intrinsic times of order 10'^ years, too short to ex- 
pect a significant evolution in the precession velocity of the 
jet). Hence a range of values of p is considered, and its ac- 
tual value is determined from the coincidence on these two 
scales of calculated values of the rate of precession of the jet, 
which exercise also determines the rate of precession itself. 

The rate of precession is evaluated as that which would 
have been observed in the absence of the lens, rather than 
that intrinsic to the source, since the component of the bulk 
velocity of the jet along the line of sight cannot be deter- 
mined (but is expected to be close to the velocity of light). 
It is calculated as described in Sections 13.11 and 13.41 over a 
range of values p = p/q. Recalling that the actual phases of 
ejection are only determined modulo mr, there is a degener- 
acy in determining the difference of phases V'2 — ipi for both 
the milliarcsecond and the microarcsecond analyses. This 
means that each separation ipi — ipj actually corresponds 
to four possible branches. Each of these branches is plot- 
ted against p in Figure |1| At this point, it is apparent that 
a value of p of ~ 0.3, and a precession rate of 0.016 radi- 
ans/day (corresponding to an observed period of 1.08 years) 
permit a coincidence of the milliarcsecond and microarcsec- 
ond scale results, so these values are adopted. To obtain the 
intrinsic source rotation period, the observed period must 
be corrected for the cosmological time dilation on account 
of the source being at a redshift of 2.5, and for the change 
in time interval on account of the jet velocity in the source. 
The first correction reduces the period by a factor of 3.5, 
and the second increases it by a factor of 1/(1 — Pios), where 
Pios = Vx/c, with c as the velocity of light. For small jet 
opening angles, f3ios ~ Vj^t/c. Assuming a value of (3 — 0.99, 
the intrinsic period is calculated to be 30.8 years. 

Deriving a period of 1.08 years from the observations 
raises a question about the authenticity of the results ob- 
tained in this paper. We have checked carefully for the pos- 
sibility of effects related to the annual motion of the Earth 
around the Sun. The paper of Jin et al. (2003), from which 
the microarcsecond scale observations are drawn, included a 
correction for the differential annual aberration between the 
NE and SW groups of images (separated by ~ 1 arcsec), 
which is typically of the order of tens of microarcseconds. 
Relative annual parallax between the NE and SW core im- 
ages is of order lO"'' /xas, too small to produce a spurious 
effect in the observations. (As an aside, another possible 
source of lensing capable of causing a differential shift in the 
location of the NE core image relative to the SW one is the 
Sun itself, which was always more than 15 degrees away from 
PKS 1830-2 11 in the observations of Jin et al. (2003), thus 
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producing a differential shift in the NE core image relative 
to the SW core image of much less than a microarcsecond) . 
Another reason that we do not consider our present result 
to be a consequence of the Earth's motion around the Sun is 
the coincidence of this result as obtained from the milliarc- 
second scale data of Garrett et al. (1998) with that obtained 
from the microarcsecond scale data of Jin et al. (2003), as 
shown in Figure 4. It is difficult to see how the annual mo- 
tion of the Earth could arrange for this, especially as the 
milliarcsecond scale data is just one actual epoch of obser- 
vation (which effectively includes two epochs for the source 
only because of the existence of two groups of lensed images 
mapping back to source events separated in time by the time 
delay of 26 days) . 

The error in evaluating the precession period is esti- 
mated via a Monte Carlo exercise that generates, for each 
of 60,000 experiments, a random selection of all the obser- 
vational errors, taking into account the cross-correlations 
between elements of the SW2NE image-to-image transfor- 
mation matrix, and then computes the desired precession 
rate for a value of p = p/q of 0.3. For the microarcsec- 
ond scale analysis, the precession rate (in radians per day 
with 1 - cr errors) is 0.0155 (-0.0166,0.0415), and for the 
milliarcsecond scale analysis, the corresponding values are 
0.0164 (-0.575,0.987). The large errors on the milharcsec- 
ond scale result are mainly due to the form of the expression 
used to evaluate the times of ejection for the two plasmons 
(see Ean.TO. 

The foregoing results are derived, as mentioned earlier, 
for an assumed perfect alignment of the line of sight with 
the jet axis. Typically, this is not the case, although it is 
expected that the alignment of the two is to within a few 
degrees for blazar sources. The effect of considering an in- 
clination of the jet axis to the line of sight can be estimated 
from expressions (1), (2) and (3) for the x, y and z veloc- 
ity components of a plasmon in the precessing jet model of 
Gower et al. (1982), where axes z and y are in the plane of 
the sky and centred on the central engine, and axis x is along 
the line of sight. In general, an additional rotation through 
an angle x of Gower et al.'s y and z axes about the x axis is 
necessary to bring the model into alignment with the right 
ascension and declination axes. Let the jet half-opening an- 
gle be (j), and let, as in earlier expressions, the angle through 
which the jet has precessed at the time of observation, as 
measured from a reference time tref, be tp. Then, with i as 
the inclination angle of the jet axis to the line of sight, the 
inferred value of vs/va from, e.g., Equation 1391 is related to 
the precession angle tp (for each plasmon considered) by: 

Sinx + {^)Cosx _ Sin(j)Sintl; 

Cosx-{^)Sinx^ Sin{i)Cos(l)-Cos{i)Sin(l)Cosi, ^ ' 

From the milliarcsecond scale data, for example, with a 
jet inclination angle of 0°, for feature 5, tp — 200°, and 
tej — tref = 125.2 days, and for feature 6, the correspond- 
ing values are —69° and 30.9 days, yielding, as mentioned 
above, a precession period intrinsic to the source of 30.8 
years. For non-zero values of i and 0, in order to obtain 
the precession rate or period in the source rest frame, the 
timescales tej — tre/ must be in each case corrected by a 
factor of 1/(1 — Vx/c), where Vx is the line of sight velocity: 

— = /? {Sincj>Sin{i)Cosil^ + Cos(l)Cos{i)) (41) 



d 













\ 1 \ 















-1 



p/q value 
Figure 4. 

Solution for tlie combination of tlie Icnsing quantity p and tfie 
rate of precession of the blazar jet (dots are used for tfie mi- 
croarcsecond scale solution, and dashed lines for the milliarcsec- 
ond scale solution). All four possible branches are evaluated for 
each solution. The coincidence of the milliarcsecond scale and the 
microarcsecond scale solution is the desired result {p/q=0.3, pre- 
cession rate = 0.016 rad/day). For p/q < —1.7 and p/q > 0.7, the 
NE image becomes negative parity, which is known to be incor- 
rect. Positive values of the precession rate are for anticlockwise 
rotation in the source (as also for the NE image; parity differences 
cause an opposite rotation in the SW image). 

With i = 8°, = 5°, X = 60° and /3 = 0.99, the precession 
period in the source frame is 21 years. A choice of i = 4.°5 
and = 3° with x and /3 as before, yields a precession period 
of about 25 years. 

3.6 Single plasmon events on the microarcsecond 
scale 

From Table 1 of Jin et al. (2003), the ratio of the NE to 
SW core image ffux densities is seen to vary as a function of 
observational epoch (c/. also the crosses plotted in Fig. |^. 
The variation shows two apparent maxima, indicating source 
variability (although variations on the timescales of a couple 
of weeks are possible with a combination of an apparently 
superluminal source and microlensing by stars in the lens 
galaxy — cf. Gopal-Krishna & Subramanian 1991 — it will 
be post-justified that this is not the case in the present cir- 
cumstance). It is tempting to model these variations in terms 
of the emission of synchrotron self-absorbed plasmons. Fol- 
lowing Expression 11 of van der Laan (1966), a series of four 
succesively emitted plasmons is considered, expanding with 
constant velocity, with an initial spectral maximum at 43 
GHz, and electron power law index of 3/2. Since the model 
is actually underconstrained by the observations, it is only 
attempted to seek a consistency between the two. For each 
image, the emission from the series of four plasmons is com- 
puted. In the case of the SW image, this curve is scaled 
by the NE/SW image fiux density ratio | k \— 1.161, de- 
rived from Section 13.31 and delayed in time relative to the 
NE curve by the time delay of 26 days. The middle two 
plasmons are constrained to have an observed separation, in 
terms of time of ejection, of (125.2 — 30.9) — 94.3 days, as 
inferred from Section f3.4l It is found necessary to introduce 
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Figure 5. 

(Solid curve): Simulated variation of the ratio of the NE image's 
flux density to that in the SW image at 43 G^z. (Crosses): 8 
epochs of data, with errors within the extent of the symbols used, 
from Jin et al. 200'i).( Dashed line): Simulated underlying flux 
density variation in the NE image core, and (dash-dot line): cor- 
responding time-delayed variation in the SW image core. Epochs 
2 and 5 of the data correspond to the commencement of the ob- 
served microarcsecond-scale 'events' discussed in Section l3.ll 

a constant core flux density, of value comparable to that of 
the plasmons, to simulate the data. The result of this exer- 
cise is shown in Figure |S] Note that corresponding to Epoch 
2, the second plasmon has just begun to initiate a rise in 
the SW curve, and again at Epoch 5, the same phenomenon 
is observed. Hence the discontinuities in Figure Q are seen 
to be a consequence of the observed appearance of a new 
plasmon in the SW image, as was stated in Section f3. II 



4 A MASSIVE BINARY BLACK HOLE 
SYSTEM FOR A CENTRAL ENGINE? 

A popular model for Active Galactic Nuclei (AGNs) with 
precessing jets is that they host at their centres massive bi- 
nary black hole systems (Begelman, Blandford & Rees 1980, 
Roos 1988, Roos, Kaastra & Hummel 1993, Villata et al. 
1998, Villata & Raiteri 1999, Romero et al. 2000, Romero, 
Fan & Nuza 2003, Rieger 2004, Maness et al. 2004). The sys- 
tem under study in the present paper lends itself to such a 
possibility. The precession period intrinsic to PKS1830-211 
has been derived from Section [3. 51 to be about 30.8 years; as 
was the case for 1928-1-738 (with a jet precession period of 
2.9 years; Roos et al. 1993) and 3C 273 (with an observed 
jet precession period of 16 years; Romero et al. 2000), this is 
too rapid to be due to geodetic precession of the jet-emitting 
black hole (Begelman, Blandford & Rees 1980), without sug- 
gesting a lifetime (~ 10^ years) for gravitational collapse of 
the binary system that is too short. One may infer, however, 
that the precession period is actually a result of the orbital 
motion of the black hole binary (as did Roos et al. 1993). In 
this case, the jet velocity is modulated by the orbital motion 
of the jet-emitting black hole. If P is the precession period, 
G the gravitational constant, and m and M the masses of 
the two black holes, then, with r as their separation from 
each other: 



For a precession period in units of 30.8 years, and scaling 
the masses by 10* Mq and the separation by 10^'' cm, 

ri6 = 6.833(P30.8)^''^(m8 -f Msf'' cm (43) 

The gravitational lifetime of the system is given by: 



tgrav = 2.9 X lO^^^^ri^j—^ yrs 



+ ^^^^ 

In this case, using Equation 1431 the lifetime of the system 
until gravitational collapse is: 



t'orav — 6.322 X 10 



8 (Mg-f 7718)1/3 



yrs 



(45) 



The motion of the core of the jet-emitting black hole, 
orbiting around its companion, would be on the scale of 
microarcseconds, undetectable in the 43 GHz observations 
used in this work. 



5 DISCUSSION AND CONCLUSIONS 

The observations of Jin et al.(2003) permit an analysis of 
the blazar source of PKS 1830-211 under atypical circum- 
stances. Firstly, changes in the core of the source on scales of 
tens of microarcseconds could be probed because of the fact 
that the source is doubly imaged within a single VLBI field 
of view, and hence uncertainties in the separation of the two 
core images could be restricted to essentially thermal noise 
limits. Secondly, there is some degree of magnification due 
to lensing, which, while its numerical value cannot be de- 
termined without the existence of a standard candle/ruler 
within the source, still provides a closer look than might 
otherwise have been possible (a typical linear magnifica- 
tion factor would be about 3). Thirdly, this lensed system 
permits a detailed study of a very distant blazar (redshift 
Zs = 2.51). These observations, coupled with data from ob- 
served changes in the milliarcsecond scale of structure, make 
it possible to estimate the precession period of the relativis- 
tic jet in the blazar source in PKS1830-211. The observed 
period turns out to be 1.08 years, assuming the jet preces- 
sion axis to be perfectly aligned with the line of sight to 
the source. This is one of the highest redshifts for which a 
blazar jet's precession period has been actually measured. 
For PKS1830-211, the typically assumed mass of ~ W^Mq 
or less for each element of the central powering black hole 
system appears to pose no particular evolutionary problems 
at this redshift (c/., e.g., Yoo & Miralda-Escude 2004). 

Estimating the mass of a possible central black hole 
system from a measured SED is confounded by the fact 
that PKS1830-211 is seen through the Galactic plane and 
is subject to much obscuration in the optical. Moreover, it 
is lensed, so there remains the possibility of spectral changes 
due to different scales of the source being magnified to dif- 
ferent extents on account of microlensing by stars in the lens 
galaxy. However, a method such as measuring the timescale 
of variability of gamma ray emission from the source can be 
used. Although two timetracks, one time-delayed with re- 
spect to the other by 26 days, would be observed together, 
the timescale of variability is expected to be typically of 
the order of a day or less (Liang & Liu 2003), making it 
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Figure 6. 

Plots of the gravitational lifetime of PKS1830-211 as a function of 
primary mass in the inferred binary black hole system (ordinates, 
in units of 10* solar masses) and the mass ratio of the two black 
holes (abcissae). The plots are for gravitational lifetimes of {from 
top down) , 10*, 10^ & 10^'' years respectively. These plots 
use the derived intrinsic period of precession of 30.8 years; this 
assumes a jet bulk velocity of 0.99c. 

unnecessary to correct for the time delayed track. Then, fol- 
lowing the analysis of Liang & Liu (2003), an upper limit 
to the mass of the jet-emitting black hole can be derived. 
(Note that these authors do in fact list a black hole mass of 
IQ^--^ Mq and a calculated minimum variability timescale 
of 10* '' s in Table 1 of their paper, for the system 1830-210, 
based on an incorrect source redshift of 1 and the observed 
EGRET gamma ray flux, which would need to be corrected 
for lensing; the absolute magnification factor is an unknown. 
If it is taken as a typical value of 10, one obtains a black hole 
mass of 10* '' M0 for the source at the observed redshift of 
2.51; note also that this mass is determined under the as- 
sumption that the gamma ray flux is isotropic or unbeamed. 
This mass corresponds to a gravitational lifetime of the bi- 
nary black hole system of order 10* years (c/. Figure (SJ). 

In passing, we note that it has been assumed in the 
present work that the plasmons follow straight line paths 
along the surface of the cone formed by the jet's precessing 
nozzle (following from ballistic ejection of plasmons). On 
the milliarcsecond scale, the essentially two epochs' worth 
of data on the source represented by the single observation of 
14 July 1996 (Garrett et al. 1997), which has been used here, 
permits no investigation of a possible curvature in plasmon 
motion. Multiepoch radio VLBI studies at, say, 43 GHz, 
during an active phase for the system, should be able to 
shed some light on the question of whether we are really 
seeing ballistically ejected plasmons, or whether the heli- 
cal appearance of the jet is accompanied by motion of the 
observed features along helical paths, in which case by the 
time the jet has reached milliarcsecond scales, observational 
evidence of curvature in the motion of some features might 
be expected. More sensitive than the observed motions of 
individual features perhaps would be information from the 
observed light curves in each image. Features moving along 
helical paths would show variable relativistic beaming if the 
jet axis is even slightly inclined to the line of sight (c/., e.g., 
Schramm et al. 1993), exhibiting a quasi-periodicity in the 



Table 1. 



Data from the observations of 14 July 1996 (Garrett ct al.l998) 



Feature 


Flux Density 




Position w.r.t. 


core imago 




(mJy) 




S,A or s)(mas) 


Al (core) 


191.7 ± 1.4 




( 0.000, 


0.000) 


A2 


58.9 ± 1.5 


(- 


-0.269 ± 0.002, 


0.299 ± 0.002) 


A3 


8.9 ±1.4 


( 

\ 


0.154 ± 0.013, 


0.663 ± 0.012) 


A4 (knot) 


36.4 ± 1.4 


( 


0.980 ± 0.003, - 


-0.150 ±0.003) 


A5 


19.7 ± 1.5 


( 


1.336 ± 0.007, - 


-0.334 ± 0.007) 


A6 


11.6 ± 2.6 


(- 


-0.582 ± 0.033, - 


-1.502 ± 0.038) 


Bl (core) 


189.2 ± 1.6 




( 0.000, 


0.000) 


B2 


19.5 ± 1.5 


( 


0.279 ± 0.005, - 


-0.079 ± 0.006) 


B3 


24.1 ± 2.4 


(- 


-0.266 ± 0.012, 


1.008 ± 0.016) 


B4 (knot) 


18.7 ± 1.4 


(^ 


-0.319 ± 0.005, 


0.655 ± 0.005) 


B5 


43.6 ±3.1 


(- 


-0.536 ± 0.015, 


0.352 ± 0.009) 


B6 


16.4 ±1.5 


( 


-0.159 ± 0.007, - 


-0.344 ± 0.007) 



intensity of emission. (Note, however, that the strict linear- 
ity apparent in Figure^spanning three and four epochs is a 
more reliable indication that on the scale of tens of microarc- 
seconds at least, the motion of the plasmons in the present 
study show linear velocities, with little sign of curvature). 

Lastly, there is a question regarding the plasmon life- 
times. Are the plasmons we 'see' in Jin et al. (2003) on the 
scale of tens of microarcseconds likely to be the progenitors 
of features such as have been seen on the scale of milliarc- 
seconds in Garrett et al. (1998)? We think this is unlikely; 
longer-lived plasmons would not permit a simple analysis 
such as has been possible in this work, with almost single 
plasmon events to be observed on the scale of tens of mi- 
croarcseconds. It would appear that the plasmon emission 
mechanism is of a highly variable nature, only occasionally 
producing long-lived plasmons such as were seen in Garrett 
et al. (1998), but for most part producing relatively short- 
lived plasmons in a quiescent phase. 
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7 APPENDIX: THE MIXED ADJUSTMENT 
MODEL 

From Section 13.21 the five constraint equations F, Equa- 
tions 1181 1191 and the three implicit in equation 1211 link 
the four unknowns, the transformation matrix elements Ti 
{i — 1...4) (constituting vector X in the present analy- 
sis), with sixteen observational quantities (positional separa- 
tions, from their respective core images, of the eight features 
A4(B4), A2(B2), A5(B5) and A6(B6), as seen in both the 
NE and the SW images. These sixteen quantities constitute 
the vector L in this section. The positions of the observa- 
tional features, with errors of fit, are listed in Table [T] The 
five equations F overdetermine the quantities in X, hence 
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a numerical method, the Mixed Adjustment Model, is em- 
ployed to solve for X. Note that the observations and the 
unknowns are implicitly related. The treatment here follows 
that described in A.Leick (1995), Section 4.4. 

Let in the course of the application of this method, the 
16 adjusted observations be denoted by La, and the adjusted 
vector of unknowns be denoted by Xa. The mathematical 
model is then: 



F(L«,Xa 



: 



(46) 



Let Xo be a set of known but approximate values of X, and 
denote the vector of observations by L(,. Then Equation 1461 
may be written as; 



F(Lb + V,Xo + X) =0, 



(47) 



where V and X are the residuals at each stage of iteration. 
This essentially nonlinear form is linearized about the point 
(L(,, Xo): 



5B16 1 
where: 



eVi + 5A4 4X1 + 5W1 



9F 



0, 



" 9L ' f'^O'^t' 
5W1 =F(Lb,Xo) 



(48) 

(49) 

(50) 
(51) 



Now a least-squares estimate of X is sought, with the tar- 
get function V"^PV, where P is the 16 x 16 weight matrix 
constructed from the covariance matrix of the observations, 
o-qSl^, where is the a priori variance of unit weight. The 
target function is minimized subject to the constraints im- 
posed by the linearized mathematical model. Equation 1481 
with a vector of Lagrangian multipliers, K. Hence, one min- 
imizes with respect to V, K and X, the function: 

(^(V,K,X) = V^PV-2K^(BV + AX + W), (52) 

and the three constraint equations arising from this process 
are used to calculate V,K and X. Since the actual math- 
ematical model is non-linear, it is necessary to iterate to 
a solution the actual values of V, K and X. The iterative 
process is said to converge if: 



(V'^PV). - (V^PV)(,_i) |< e. 



(53) 



where e is a small positive number (here taken as 10~*). In 
the present case, the calculation converged in five iterations, 
with a final value of V'^PV=9. The estimated a posteriori 
variance of unit weight, a^, is therefore V^PV/(5-4)=9. 
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